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We study the potential impact of improved future supernovae data on our understanding of the 
dark energy problem. We carefully examine the relative utility of different fitting functions that can 
be used to parameterize the dark energy models, and provide concrete reasons why a particular choice 
(based on a parameterization of the equation of state) is better in almost all cases. We discuss the 
details of a representative sample of dark energy models and show how future supernova observations 
could distinguish among these. As a specific example, we consider the proposed "SNAP" satellite 
which is planned to observe around 2000 supernovae. We show how a SNAP-class data set taken 
alone would be a powerful discriminator among a family of models that would be approximated 
fSJ ' by a constant equation of state for the most recent epoch of cosmic expansion. We show how this 

, family includes most of the dark energy models proposed so far. We then show how an independent 

■ measurement of Sim can allow SNAP to probe the evolution of the equation of state as well, allowing 

(N ■ further discrimination among a larger class of proposed dark energy models. We study the impact 

$— ( , of the satellite design parameters on this method to distinguish the models and compare SNAP to 

alternative measurements. We establish that if we exploit the full precision of SNAP it provides a 
I very powerful probe. 
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I. INTRODUCTION 
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One of the most challenging problems in modern cosmology is to provide an explanation for the recently observed 
accelerated expansion of the universe . These observations have reopened the quest for the cosmological constant 
which was introduced by Einstein but later abandoned [|j and infamously cited as his greatest blunder ||]. The 
cosmological constant can be considered as new kind of "world matter" ||^ and be identified with the energy density of 
' the vacuum explaining and computing it in terms of particle physics has been largely unsuccessful [^,0 because 
it is very difficult to explain the small vacuum energy density of 10~^^''Mp[ within fundamental physics; typically it 
Q ' is either much larger or exactly zero. 

^ I In recent years, the type la supernovae (SNe) as standard candles have been used to measure the distance-redshift 
' relation in the universe providing evidence for an energy component in the universe which behaves like a cosmological 



constant This means the pressure of this component is negative and it appears to be dark in a sense that it 

is not recognizable by direct observation The Supernovae Cosmology Project (SCP) [^]J^ found evidence for a 
positive cosmological constant on the 99% level. These findings seem to be confirmed if one combines the most recent 
cosmic microwave background (CMB) radiation data from the BOOMERanG (Balloon Observations of Millimetric 



Extragalactic Radiation and Geomagnetics) [y_2-15|, MAXIMA (Millimeter Anisotropy experiment IMaging Array) 
p6|-p8| and DASI (Degree Angular Scale Interferometer) [^,^ experiments with observations of rich clusters pl| , p2| . 

With these observations we need a deeper understanding of the cosmological constant and attempts have been made 
to explain the missing energy as the energy density in a scalar field, which only interacts with the other fields via 
gravity. This field is rolling slowly down a potential or gets trapped in a local minimum ||23|-p0| . Therefore the vacuum 
energy of the universe becomes important for its evolution and the expansion begins to accelerate, generalizing the 
concept of the cosmological constant. Attempts have been made to connect this field to fundamental physics [ p9|38 



and resolve the problem of fine tuning of initial conditions p6| . The problem is that there is a plethora of models 
which can describe the observed expansion, but with the current available data it is not possible to distinguish between 
most of them. 

To improve the observational situation a satellite mission - the "SuperNovae Acceleration Probe" (SNAP) - Q 
and other dedicated SNe surveys have been proposed . This satellite may observe about 2000 SNe within two years 
and therefore increase the number of SNe by a factor of 25. In this paper we present the details of a representative 
sample of dark energy models, discuss how the use of SNe as standard candles can distinguish the different models 
and what can be established about the equation of state of the dark energy component. Current upper bounds from 
SNe observations on the equation of state are < —0.6 | pT| , ^j45| . In order to reconstruct not only the constant 
contribution to the equation of state it is convenient to fit the SNe magnitude - redshift relation with a continuous 
hmction @-g|,|6|-|8|. 
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The main purpose of this paper is twofold: First we compare the quahty of two different fits to the luminosity 
distance - redshift relation, where we emphasize the importance of a "good" fit in order to draw conclusion about 
the quality of the experiment. Secondly we analyze how SNe observations can constrain the equation of state factor 
and content of dark energy in the universe. A whole section presents the different dark energy models which we use 
for this analysis. 

The outline of the paper is as follows. In section |l| we describe the current situation of the SNe observations, 
in section |l| we describe briefly the specifics of the proposed SNAP satellite mission, in section IV we introduce 
the commonly studied dark energy models, the parameters we choose for them and their cosmological evolution. In 
section ^ we discuss how to reconstruct the equation of state by expanding the equation of state factor as a power 
series in redshift and fit for the expansion coefficients. In section VI we discuss the impact of the experimental design 
and prior constrai nts on the matter content. Finally, section VIl presents alternative measurements, before we draw 
our conclusions in VIII. 



II. THE CURRENT OBSERVATIONAL SITUATION 

In this paper we concentrate on the results of the distant type la SNe observations [0-^ and mention other 
indications just briefly. The SCP and the High-Z Search Team used bright type la SNe as standard candles. These 
objects are thought to be thermonuclear explosions of carbon-oxide white dwarfs |50|-[52[. The correlation between the 



peak luminosity and the decline rate of the luminosity of the SNe 1 53 - 5^] makes it possible to estimate its magnitude 
and with spectral information about the host one can determine its redshift. The correlation can be quantified by 
the drop in magnitude 15 days after the peak luminosity is reached. The SNe observations by the SCP are calibrated 
using the "low" redshift Calan/Tololo survey which revealed that they have an excellent distance precision of 
Cmag = 0.15 mag for the magnitude and therefore can, in fact, be used as standard candles. 
The apparent bolometric magnitude is given by 

to(z) = Af + 51ogdL(z) + 25, (1) 

with M the absolute bolometric magnitude and dh the luminosity distance, which is usually defined with distances 
in units of lOpc. However, cosmological distances are measured in Mpc and therefore there is an additional term 
5 log 10^ = 25 in eqn. (|l|). Furthermore, the luminosity distance depends on the cosmological evolution and hence 
on the cosmological parameters and is defined by c?l = C/AnJ-, where J- is the measured flux and C the absolute 
luminosity of the object. The luminosity distance can be expressed in terms of the coordinate distance di,{z) — 
(1 + z)r(z). As mentioned before the CMB data of BOOMERanG, MAXIMA and DASI in combination with 

2dF observations indicates strongly that the universe has a flat topology and therefore we concentrate here on 
flat cosmologies where the coordinate distance is given by 



In expression (Q) the quantities m, M and c^l depend on the Hubble parameter Hq. From eqn. (H) we see that 
dL ~ so we can rewrite m{z) in the following way 

m(z) = M + 51ogDL-51ogi?o + 25, (3) 

where we have defined Dl = H^di^. For low redshift SNe we can use the linear Hubble relation 

m(z) = M + 5 log cz - 5 log Hq + 25 

= + 5 log cz , (4) 

where we have defined the magnitude "zero point" Ai = M — 51ogi/o + 25. Theoretically, this quantity can be 
determined by the survey, but in practise this is just a statistical nuisance parameter which is marginalized to 
estimate the cosmological parameters, so that Di^{z) can be estimated without explicit knowledge of Hq fl . In fig.|l| 
we plot the effective bolometric magnitude mf data points of the SCP and Calan/Tololo survey pj\ as well as 



the curves m{z) from the theoretical models we study in section IV. The effective magnitude refers to the apparent 
bolometric magnitude which has been corrected by the lightcurve width-luminosity correction, galactic extinction and 
the K-correction from the differences of the R- and B-band filter If we just allow a cosmological constant, dark 
and baryonic matter content in the universe and assume a flat cosmology, which seems to be confirmed by recent 
CMB and large scale structure observations p^ , p^ , p7|j58| , the best fit values are roughly fim = 0.28lo;Qg and therefore 
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flA = 0.72 §. This is in agreement with the analysis of the High-Z Search Team of = 0.24 and flj^ = 0.76 
The low matter density is confirmed by several observations; by the evolution of the number density of rich 
clusters mass estimates of galaxy clusters, either by the Sunyaev-Zel'dovich effect or through measurements 
of the X-ray flux pi| and also the shape of the matter power spectrum |plt|. Visually, one can hardly distinguish 




FIG. 1. The Calan Tololo (open circles) and SCP data points (solid circles). The curves correspond to the theoretical models 
discussed in section 

the different models in fig.|l| in a sense that they all seem to fit well. In order to be able to tell more we plot the 
magnitude difference Am(z) = m(z) — m/^{z), where niA is the apparent magnitude of a fiducial model with = 0.7 
and fini = 0.3. We show the relative magnitude in fig.|[ We can clearly distinguish most of the models of section IV , 
however we already realize the problem with the current observational situation that the available data can not really 
differentiate between the particular models. If we want to distinguish the models in the future we have to be able to 
achieve much smaller errorbars than the SCP data has. In the next section we describe the SNAP satellite project 
which is able to achieve this goal. 



III. THE SUPERNOVAE ACCELERATION PROBE - SNAP 



In order to improve the current observational situation significantly a new satellite, the SuperNovae Acceleration 
Probe - SNAP has been proposed which will be dedicated to the observation of SNe. The SNAP satellite is 
equipped with a 2 meter telescope with a optical imager, a ID' near-IR imager, and a three-channel near-UV-to- 
near-IR spectrograph. Every SNe at z < 1.2 will be followed as its brightens and fades. The wide-field imager makes it 
possible to find and follow approximately 2000 SNe la in two years, and the 1.8— to 2.0— aperture of the mirror allows 
this data set to extend to redshift z = 1.7. Furthermore, systematic uncertainties will improve considerably compared 
to the current situation. The uncertainty due to the Malmquist bias, the fact that the most distant SNe are only 
the ones with large intrinsic brightness and therefore represent a very biased sample in brightness, will also improve 
since each SNe will be observed 3.8 magnitudes below peak brightness. For a large subsample spectral time series 
and cross-wavelength flux calibration will reduce the uncertainties from the K-correction and cross-filter calibration. 
The Sloan Digital Sky Survey (SDSS) the Space InfraRed Telescope Facility (SIRTF) || observation and SNAP 
spectra of host galaxy subdwarfs will improve the systematic uncertainty due to the Milky Way Galaxy extinction. 
The uncertainty due to gravitational lensing by clumped masses will be averaged out due to the large statistics. The 
error due to extinction due to ordinary dust outside the Milky Way will be reduced due to the cross-wavelength 
calibrated spectra. Also the uncertainties due to Non-SNe contaminations will decrease. Gray dust uncertainties 
can be addressed due to large observed redshifts, with z > 1.4, and with broad wavelength measurements into the 
near-IR. Due to the large sample size and detailed lightcurve and spectral information, SNAP will provide sufficient 
data to measure second order effects like uncorrected evolution of the SNe. These systematic errors lead to an absolute 
uncertainty of CTsys = 0.02 mag at redshift z = 1.5, while the statistical calibrated uncertainty is Umag — 0.15 mag 
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which corresponds to approximately 7% uncertainty in the luminosity distance. The redshift coverage of the SNAP 
satellite within two years is shown in table ^ The numbers in table | are based on the observed rates of SNe out to 



redshift interval 


2 = - 0.2 


z = 0.2 - 1.2 


z = 1.2 - 1.4 


z = 1.4 - 1.7 


number of SNe 


50 


1800 


50 


15 



TABLE I. SNAP specifications for a two year period of observations, with a statistical uncertainty of (Jmag = 0.15 mag and 
an uncertainty limit of asys = 0.02 mag at redshift z = 1.5. 



redshift z = 1.7. In the redshift interval z — 0.2 to z = 1.2 we assume that SNAP will observe a 20 field within 
two years. At high redshifts there are many more SNe, but SNAP will not have the time for a spectroscopic follow 
up on all of them. Likewise, at the lowest redshifts there will be more type la SNe, but the limiting factor here is 
the sky coverage of SNAP. We use these specifications and simulate the SNAP experiment assuming a background 
cosmology with fim = 0.3 and J^a — 0.7. In the simulation we assume a Gaussian distribution of the uncertainties 
and an equidistant sampling of the redshifts in the four ranges. We further neglected the errors in redshift, since 
they are expected to be of the order Sz — 0.002 and therefore relatively small. In fig.^ we show the results of this 
simulation. For plotting purposes we bin the data points so the resulting uncertainty is a = 0.02. The number of data 
points in one bin, A^bin, is given by A'bin < "'mag/^'sys- However, the realistic situation is a bit more tricky since the 
systematic error is drifting from cTgys = at z = to cTsys = 0.02 at 2 = 1.5. For the discussion of systematic errors 
we assume a linear drift and use CTsys = -ZCTsys/l-S. In the right plot in fig.|^ we also bin the SCP and Calan/Tololo 
data. One clearly recognizes in fig.0, that one can distinguish some of the models with a SNAP type observation. 




FIG. 2. The relative magnitude with respect to a cosmology with fim = 0.3 and Qa ~ 0.7. 
simulated with this cosmology. The solid triangles are the btnned data points with errorbars from 
as in table |. We have not plotted the data in the redshift interval 2 = — 0.2 for the SNAP 
Calan/Tololo (open circles) and SCP data points (solid circles) are not binned and in the right 
correspond to the theoretical models discussed in section ^ and the key to this curves is the 
dot-short dash line is a cosmological constant model with 57a = 0.6 and the thick short dash - 
r^A = 0.8. The thick long dash line is the "Standard Cold Dark Matter" model with f2m = 1.0 
the current data. 



The SNAP data points are 
the SNAP type specifications 
experiment. On the left the 
figure they are. The curves 
same as in fig.^ The thick 
long dash line a model with 
which is clearly ruled out by 



while the current data does not allow any differentiation. In the next section we will present the dark energy models 
we studied and then in section ^ we will quantify how SNe observations can distinguish these models. 



IV. DARK ENERGY MODELS 



As mentioned in the introduction a possibility to generalize the concept of a cosmological constant is by introducing 
a scalar field which only gravitationally interacts with the other fields. The dark energy field is supposed to slowly 
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roll down the potential or is trapped in a local minimum. This leads to a vacuum dominated state of the universe 
which hence leads to an accelerated expansion. The energy density of the field is given by its kinetic and potential 
component, 



1 -a 



while the pressure is given by the difference, 



P0 = or + v{c^) , (5) 



P0 = - y{<f) ■ (6) 



2 

Note that we assume that the field is homogeneous on large scales. The proportionality factor 

w^ = ^ , (7) 
P4> 

in the equation of state, = Wfj^pcj^, \s — —1 ii the kinetic term (^^/2 is negligible. This is exactly the equation 
of state for a cosmological constant term. In this paper we study the behaviour of the magnitude - redshift relation 
and therefore, we have to solve the Friedmann equation 



a\ 1 



a I 3 



Pother + + 



(8) 



where we have used the Planck mass Mpi = 2.44 x lO^^GeV as a unit, a is the scale factor of the Robertson- Walker 
metric and Pothor is the total energy density of the other contributing fields or energy components, like dark and 
baryonic matter and radiation. The evolution of the dark energy field is given by the field equation 

4> + iH4> + V'{(t>)^Q, (9) 

with V {(j)) — dV /d(j). If V{4>) is approximately constant and the other energy components are negligible, the solution 
for the scale factor is a ^ exp[A/y t] and hence, the expansion of the universe is accelerating. This is the same concept 
as inflation ||6^-|66||, which also exploits the rapid expansion rate. However, in the context of the cosmological constant 
and dark energy we are interested in solution where the universe is vacuum dominated only in recent times and not 
in the early universe as in inflationary models. 

There are two possibilities to neglect the kinetic energy 0^/2: either the field rolls down very slow the potential 
"hill" or it is trapped in a local minimum, which is illustrated in fig.^ On the left side we plot an exponential potential 
as discussed in [p5]-|2^ 
model proposed in 



which gives rise to a slow roll of the dark energy field. On the right side we plot the 
lich is a potential with a local minimum. In fig.^ we show the evolution of the densities 
relative to the critical density pc = l/3H^{z) for the model described in |^0|. This example shows the generic feature 
that the universe first is radiation dominated (dashed line), then matter dominated (dotted line) and finally becomes 
dominated by dark energy just before the present day (solid line). 

In the following discussion we present the dark energy models we will use to test the opportunities of future SNe 
observations. The specified parameters of the models below all lead to Hq = 65 km/sec/Mpc, flm = 0.3 and fi^ = 0.7, 
where is the relative energy density of the dark energy component today. 

Pure exponential - This potential appears naturally in compactified higher dimensional Kaluza-Klein theories as well 




as in certain Supergravity models. This model was discussed in the context of a dark energy field in p5|- p8| , plf ^ 73 
and is given by the potential 

y(</)) = Kje-^^ (10) 

For a range of parameters and initial conditions this solution exhibits an attractive behaviour whereby the field 
tracks the dominant component of the background cosmology, i. e. matter or radiation. However, in order to fulfill 
the observational constraint of the SNe experiments the model parameters have to be chosen from the transient, 
non-attractive, branch and the model needs fine tuning of initial conditions: in our discussion we used the values 
Vo = 10~'^^"M^i, A = IMp/, 0(0) = 0.135Mpi and 0(0) = 0. The magnitude - redshift relation for this model 
corresponds to the thick long-dashed line in figs. Note that this model is slightly below the zero line of the A 

model in fig.|| and just outside the Icr errorbars of the SNAP satellite. 

Pseudo Namhu-Goldstone Boson (PNGB) - This potential can arise as potential energy of very light axions if the 
U{1) Pecci-Quinn symmetry is broken p9| , |30| , |3^ , p9| . The potential is given by 
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FIG. 4. The evolution of the densities relative to the critical density for the trapped minimum model |40||. The long dashed 
line is Qi, the energy density in the radiation, the dotted line f2m, for the matter fields, and the solid line 57^, the dark energy 
contribution. 



-M^ [cos ((/.//) + 1] , (11) 

where M is of the order of a very Ught neutrino mass {M ^ 0.001 — O.OleV) and / is the symmetry breaking scale 
(/ ~ 10^'°^ — lO^^GeV) . For a wide range of parameters the model can behave hke a pure cosmological constant and does 
not require fine tuning. However, we studied a parameter branch where the equation of state factor oscillates, as seen 
in fig. ^ We use this particular setup, because it will be used later to illustrate a case in which the reconstruction of the 
equation of state could be troublesome. We used M'^ = 1.001 x 10~^^°Mp[ and / = O.lMpj. In this parameter branch 
it is also necessary to tune the initial conditions to fulfill the observational constraints to ^(0) — 1.184 x 10~'*Mpi and 
4>{0) = 0. In the magnitude - redshift relation of figs .|l|,||j5| the thick solid line represents to this model. We note the 
oscillatory nature of the potential is also observed in the apparent magnitude m(z). For this choice of parameters the 
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model is almost ruled out already by the current SCP and Calan/Tololo data. 

Cosmological tracker solutions - These solutions are a generalization of the attractor behaviour of the pure exponential 
potential pTj . The potentials have a functional form and the most studied examples are the inverse tracker 

potential 

^W = ^' (12) 

and the exponential tracker potential 

= Af^e*^^/'^ . (13) 

The notion of tracker solutions refers to the fact that these solutions evolve on a common evolutionary track inde- 
pendent of the initial conditions |^,^. The inverse tracking potential is motivated by supersymmetric QCD. The 
common feature of these models is that the density in the dark energy field at late times dominates over all the 
other energy contributions and therefore the expansion of the universe begins to accelerate. The cosmic coincidence 



problem |35| is the fact that one still has to adjust the parameters of the model to determine the time when the 
dark energy component begins to dominate. However, as mentioned above, the initial conditions are almost arbitrary. 
For the inverse tracker potential we used the parameters M = 2.11 x lO^^^Mpi and a = 6. In figs .|l|,||,H the inverse 
tracker model corresponds to the thick dotted line. This model seems also to be marginally disfavoured by current 
data as evident from the right panel of fig. |[ The only parameter to adjust for the exponential tracker potential is 
M — 9.09 X 10^'^^Afpi and this model is plotted as a thin long dashed line in figs.0,||j| and behaves in the shown 
redshift range almost entirely like the cosmological constant model. 

Supergravity potential - This model is inspired by supersymmetry breaking in type I string theory and Supergravity 
p7| , p8| with the potential given by 



1 / 



2 \Mpi 



(14) 



Since supersymmetry breaking should occur above the electroweak scale and in order to avoid fine tuning of initial 
conditions, the parameters of this model have to fulfill the constraints a > 11 and M > lO^^Afpi. These requirements 
seem to lead to a "unnatural" way of supersymmetry breaking | |67t but nevertheless this model is rare in that it is at 
least related to a fundamental theory, and recent work shows that SUGRA may prevent this type of difficulty [ p8[ . 
For small values of cj) the exponential in eqn.(^) is approximately constant so at early times the evolution behaves 
like an inverse tracker model and has, therefore, all the advantages of the tracking solutions. The parameters we chose 
for our discussion are AI — 1.611 x 10~^A'/pi and a = 11 and the model is plotted as a thin solid line in figs.|l|,||,|[ 
We recognize that although this model is clearly different from a cosmological constant, we can not distinguish this 
model from a pure cosmological constant with the current data. 

Exponential with polynomial and rational pref actor - The problem with the models discussed so far is that the involved 
mass scales seem not to be natural in terms of Planck scale physics. The model proposed by Albrecht and Skordis 
pO| addresses this issue by multiplying an exponential potential by a polynomial prefactor 

V{(p) = yp(0)e-^^ , (15) 

where X^((/') is chosen to be 

ypi^) = A+''-^-^- (16) 

In the right panel of fig.^we show this potential for a — 2. In this example the field gets trapped in the local minimum 
of the potential independent of the initial conditions, so no fine tuningof them is necessary. In p9| it is shown that 
this false vacuum state of the field is stable to quantum decay, while |7y] discusses the possibilities of a classical roll 
through this potential. The parameters used here are A — 8Mpi^, B — 33.989Mpi and A = O.OlMpj. This model is 
shown as a thin short dashed line in figsj^,^,^ and we recognize that it is completely indistinguishable from a pure 
cosmological constant. It is possible to generalize the polynomial prefactor and allow rational functions and there 
might be a possibility to connect the dark energy potential to the interaction of two separated 3-dimensional branes 
from string theory |7l],p9| . A promising candidate for such a potential is 



7 



(17) 



where S regularises the singularity at (j) = B. As with the trapped minimum model this potential also has the 
feature that the field gets trapped in a false vacuum state. We take the parameters of the model to be A = 8Mp^, 
= B — 35.1628Mpi and 5 = O.OlMpj. This brane model is depicted as a thin short-dash dotted line in figs .|l|,H,|| and is 
also not distinguishable from a pure cosmological constant model. A model with different parameters but very similar 
behaviour is discussed in El|. In both models the involved parameters are of order 0(1) in units of the Planck mass 
Mpi. The parameter B needs to be adjusted that the field gets trapped in the local minimmn at the right time to 
account for the observed density in the dark energy field of = 0.7 today. 

Two exponentials - This type of potential could arise in string theory as a possible result of Kaluza-Klein type 
compactification and is given by 

V{(l)) = Vo [e^^ + (18) 

and there is no fine tuning problem of the initial condition s |7^ . The parameters in this model are chosen to be 
Vo = 8.2 X lO^^^^Mpj, A — 20Mpi^ and /3 — 0.5Mpi. In |72] other possible parameter choices are discussed. In 
figs.|l|,H,^ the model is drawn as a thin short-dashed - long-dashed line and we note that this model's apparent 
magnitude evolves almost like a cosmological constant in the observed redshift range. 

Periodic Potential - The common feature of most of the models discussed so far is that the parameters have to be 
adjusted in a way that the dark energy component only becomes dominant today, which means we live in a special 
epoch. The only exception is the pure exponential ]27| , |34 | which has an attractor behaviour and follows the dominant 
component of the background component, that is matter or radiation. However, the pure exponential models in the 
attractor branch of the parameter space are ruled out by Big Bang Nucleosynthesis (BBN) and the SNe observations. 
A sinusoidal modulation of the pure exponential can resolve this problem |74| and such potential is given by 

V{^) ^Vo[l + 5 sin {(3(j})] e-^^ . (19) 



There is only an adjustment of the parameters necessary to fulfill the BBN constraint and the parameters used in our 

A = 4.0M-\ S = 0.98 and (3 = O.SIM"/. 



discussion are Vq = 2.55 x lO^M^i, A = 4.0M"\ S = 0.98 and /? = 0.51M~^. This model corresponds to the thick 



short-dashed line in figs. 
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FIG. 5. The redshift evolution of equation of state factor = P<t>/p<t> for the discussed models. The thin short dash line is 
the trapped minimum model, the thin dot - short dash line is from the brane inspired potential, the thin short dash - long dash 
line from the potential which involves two exponentials, the thick short dashed line from the periodic potential, the thick long 
dashed line from the pure exponential, the thick solid line from the Pseudo Nambu-Gotu Boson potential, the thin solid line 
from the Supergravity inspired potential, the thin long dashed line from the exponential tracker solution (underneath w = — 1), 
and the thick dotted line from the inverse tracker. 
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In fig.||we show the evolution of the equation of state factor of the dark energy component. We recognize that 
most of the models have a smooth behavior, apart from the PNGB model. For this model oscillates between —1 
and 1. We have now a fairly representative sample of dark energy models though it seems impossible to include the 
rapidly increasing number of all suggested models. Two classes of models which are entirely missing in our discussion 
are the one where the dark energy field is non- minimally coupled to gravit y | |7 5[-[77| , where the field is directly coupled 
to matter and where the dark energy field is kinetically driven ||80| , |8l[] . The reconstruction of the equation of 

state for these models is discussed in In the following we will show which of the presented models can be 

distinguished by SNAP type SNe observations. 



V. A FIT DESIGNED TO RECONSTRUCT THE EQUATION OF STATE 

In order to distinguish between different dark energy models we have to be able to quantify how well SNe observations 
with a SNAP type experiment can map out the magnitude - or luminosity distance - redshift relation. There has 
been a suggestion to fit the luminosity distance by a polynomial . This fit is motivated by the need of a smooth 
function to reconstruct the equation of state factor of the dark energy component and by its obvious simplicity. More 
recent work suggest to fit the luminosity distance by a rational function with three free coefficients |^^. This fit has 
the advantage that in extreme cases it behaves like the analytical solutions for a pure cosmological constant cosmology 
or a completely matter dominated universe. There have also been other suggestions with different fitting functions 
psf and it may to be possible to fit directly for the evolution of the dark energy density 

The polynomial fit of the luminosity distance is defined by 

N 

di^iz) =Y,c,z\ (20) 

■i=0 

where we will truncate the power series at an appropriate N. Since the luminosity distance c?l(0) = 0, for all 
cosmol ogic al models, we can set cq — 0. In order to study this fit we use the proposed dark energy models from 
section IV and create data sets with the SNAP type specifications from table | with a Monte Carlo simulation. We 
assume that the errors 5m in the magnitude are Gaussian distributed with a zero mean and a variance of Umag- 
Furthermore we impose an equidistant sampling in redshift, which seems to be the optimal sampling fs^ , and neglect 
the uncertainties in redshift space. The simulation is repeated and the fitting procedure A'stat = 1000 times to obtain 
the appropriate statistics and we find that the distribution of the coefficients Ci is Gaussian. 

We will now discuss how to reconstruct the equation of state factor from the measured magnitude or luminosity 
distance. The conservation of energy in the dark energy component yields 

e±=-:^H{l + w^), (21) 

with H = a/ a — (pm + p^)/3, a^^ — 1 + z and the definition of the coordinate distance r{z) in eqn. (||) we obtain ]4^ ] 

1^ l + z 3angg(l + z)V2^ 

I + Wa, — 5 7 , 22) 

' 3 n„^H^{l + zf~^ 

where r' denotes the derivative of the coordinate distance r with respect to redshift z. Since the coordinate distance 
is r{z) = dh{z)/{l + z) we can calculate the derivatives from the fit in eqn. (^0|), although we will need to quantify 
the matter content ilm in order to do this. Since the errors in the coefficients Ci are Gaussian we can calculate the 
error in the reconstructed by ordinary Gaussian error propagation 

ij 

with aij the covariance matrix of the simulated sample of c^. 

In fig. ^ we show the reconstructed equation of state factor for different values of N. The background cosmology 
is a cosmological constant model with ilm — 0.3 and JIa — 0.7. In the top right panel we plot the mean values of the 
reconstructed w^. We note that for iV = 3 fit (dotted line and top right panel) the mean value does not represent 
the theoretical well. This is because in order to reconstruct we need the second derivative of the coordinate 
distance r{z) (eqn. E^). The coordinate distance is already reduced by one order in z compared to the luminosity 
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FIG. 6. The reconstructed equation of state factor w^f, for the cosmological constant model with the theoretical value of 
11)0 (2) = — 1. In the top left panel we show the mean values for the reconstructed with the dotted line for the A'^ = 3 
polynomial fit, the solid line the A'^ = 4 and the dashed line the A'^ = 5 fit. The top right panel is the A'^ = 3 fit with the shaded 
region representing the la uncertainty levels. The lower left panel is the same plot for A'" = 4 and the lower right panel for 
N = 5. 

distance so the iV = 3 fit might not represent sufficiently the evolution in r{z) to produce a second derivative in 
r{z) which represents the theoretical value at least roughly. The iV = 4 (solid line) and N = 5 (dashed line) fit 
however reproduce, at least in the relevant redshift range, the theoretical value of w,p to a satisfactory level. The 
= 5 fit naturally reproduces the better then the = 4 fit, however the errorbars for the iV = 5 fit (lower right 
panel) are much larger then for the TV = 4 fit (lower left panel). A fourth order polynomial, therefore, yields the best 
reconstruction for with the SNAP type specifications from table |. We note, however, from the lower left panel in 
fig. ^ that for the 4th order polynomial fit we can only reproduce in the required range if —1.3 < Wcf, < —0.7, at 
the Icr level. 

As we saw in the previous section it seems nearly impossible to make any predictions about the equation of state 
factor with the SNAP type specifications. However, the problem could be the polynomial fit and not the SNAP 
experiment. If we assume for example a constant and try to reconstruct it with a polynomial fit of the luminosity 
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distance, we can show analytically that no finite order polynomial can do this exactly. Even in the extreme case of 
just a pure cosmological constant model with J^a = 1 or the SCDM model with = 1 the polynomial can not fit it 
exactly. This problem was recognized in ^] and they suggested to use Pade approximants or even splines. In ^] 
a rational function is used which at least allows that the extreme cases of = 1 and f^m = 1 to be fitted with an 
exact relation. This method has been improved due to the introduction of a more complicated rational function with 
more free parameters [^ . A problem with all these reconstruction methods is that wc need, the matter density 
as another input parameter as seen in eqn. (p2[). A further problem is that in order to reconstruct we have to 
calculate second order derivatives in r(z) which will in general, increase the errorbars on the reconstructed equation 
of state factor. We propose a fit which has been discussed recently |8^-|8^ allowing one to read ofi^ the equation of 
state factor directly with no reconstruction as in eqn. (|2^) and also includes the possibility to fit for fl^ in a more 
natural way. 

We expand the equation of state factor into its rcdshift evolution 



N 

E 

i=0 



w., 



where we chose the expansion in (1 + z) for computational convenience. With this expansion and eqns. 
obtain the luminosity distance in a flat universe 



c{l + z) 
Ho 



: dz' . 



m.^ + n^il + zf" exp<{3E ■ 

1=1 



N 



(24) 



we 



(25) 



[1 + z' 



with 51^ = 1 — flni- We note that for Wi — for i > 1 and wq = — 1 we obtain the standard result for the cosmological 
constant. In the following discussion we take fim = 0.3 fixed and later we examine the fit with different prior 
information on Q^^. 

The fit is done by minimizing the function 



am}) = E 



fe=0 



di^jzk) - dl^Zk) 
Sdi^{zk) 



(26) 



with ddi^{zk) ~ cr,„ag'^L(-2fc) lii(10/5), the uncertainties in di^{z) as the weight on the particular data points. The 
sum runs over the whole redshift range from table |, with Nz the overall number of measurement redshifts. Again 
we assume that the redshift range is split into equidistant samples in the four ranges from table ^ We minimize 
this rather complicated expression with the MINUIT routine from the CERN program library which also delivers the 
covariance matrix on the parameters. In fig. ^ we plot the fit for different values of N for the periodic potential as 
the dark energy model. We choose this model because we see in fig. ^ that the equation of state factor for this model 
is evolving within the relevant range. For a cosmological constant model a fit with = should be the best, because 
it is exact and we do not gain any information by going to higher order. In fig. ^ we note that the N = fit, with 
~ 31, is a relatively poor fit of the theoretical values if we set a prior of flm — 0.3 (left panel). The first order 
results in a satisfactory fit with ~ 0.47 and N — 2 improves this result only slightly, with — 7.3 x 10"'^. So 
in order to study the luminosity distance c^l the first order fit seems to be sufficient. If we release the prior on 
and just constrain the matter contents of the universe by < < 1 the A^ = 1 and A = 2 fit are indistinguishable 
and fit reasonably well. It appears that the fit (eqn. p5| ) leads to more accurate results than the polynomial fit 
(eqn. |2^) with less free parameters! 

In fig. ^ we plot the relative accuracy of the different fits, for an unweighted sampling points. We recognize that 
the cubic polynomial- and quadratic ly-expansion lead to the best fit. This results holds also if we use the weights 
from the SNAP specifications. We performed this comparison for the SUGRA and the periodic potential models as 
well as for a toy model . For the cosmological constant the w-expansion results by construction to the best results 
because the fit is "exact" . We expect that this behaviour holds as well for the nearly constant models, like the two 
exponentials model. 

In fig. ^ we plot the results of the constant, A^ = 0, fit and the linear, A^ = 1, fit for a few of the sample of dark 
energy models which we have discussed (black lines) and their theoretical values (grey lines). We note that for the 
models which do not evolve much such as the inverse tracker (dotted line) and the pure exponential (long dashed line), 
the constant fit seems to be sufficient. However, for the SUGRA model (thin solid line) and the periodic potential 
(short dashed line) only the first order fit is acceptable. Note that both orders are not a good enough to fit the PNGB 
model (thick solid fine). 
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FIG. 7. The fit from eqn. (^5|) for the periodic potential as the cosmological background. We plot the fitted cZl* with respect 
to the theoretical luminosity distance di, in percent (%). The dotted line is the A'' = fit, the solid line the A'' = 1 fit and the 
dashed line the N — 2 fit. The left panel is the one with flm = 0.3 as a prior and the right one has just the constraint that 
< fim < 1. 




FIG. 8. Relative accuracy of the rational fit (dotted line), the quadratic and cubic polynomial fit (dashed line) and the linear 
and quadratic ?j;-expansion (solid line). The cosmology is taken from the periodic potential. 



In the following wc will discuss the resulting w^{z) from the fit. First, we have to know the error matrix which is 
calculated as the inverse of the second derivative of the x^-function at its minimum 



dwidw-j 



(27) 



which is a valid approximation if the function has an approximately parabolic shape around its minimum. MINUIT 
also calculates the marginalized errors on the parameters by calculating the values of the parameters for Xmin + 
with = 1. We used both methods and found that they give consistent results. The errors on are then given 
by Gaussian error propagation 
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FIG. 9. The relative magnitude plots for the A'^ = fit in the left panel and the A*' = 1 fit in the right panel. The light grey 
lines are the theoretical values and the dark lines are the fitted results. The line styles are the same as in fig. ^. 



^ dwi dwi 



(1 + ^) 



(28) 



In fig. y^we plot the resulting for different orders iV of the fit in eqn. ( pSD with the periodic potential as background 
cosmology. In the top left panel we plot the mean values and the theoretical curve (thin dotted line) . The dashed line 
is for the iV = fit and we recognize that a — wq — const, fit can not reproduce the evolving model. The = 1 
fit (solid line) already represents some evolution, so for z > 0.6 the fit becomes fairly poor. The dashed line is the 
quadratic, N = 2, fit, which leads to a better result than the linear fit. In the top right panel we show the constant, 

= fit with it's very small errorbars, in the lower left panel the iV = 1 fit and in the lower right panel the N — 2 
fit with the errorbars. We recognize that the errorbars for the iV = 1 and = 2 fit are roughly on the same level for 
z < 0.7, but then the errorbars of the = 2 fit increase rapidly. In general the errorbars on the TV = 1 fit are smaller 
then the ones on the iV = 2 fit, since there are less degrees of freedom. 

In fig. [ill we plot the periodic potential model (top left panel), the inverse tracker (top right panel) and the SUGRA 
inspired dark energy model (bottom left panel). The light shaded regions are the la confidence levels. In each plot 
we show the cosmological constant model for comparison (dark shaded region). We recognize that for the fit from 
eqn. ( p5|) we can distinguish the w^{z) evolution from a cosmological constant on the la level. We also performed the 
fits for the N — 2 approximation and got still better results as for the polynomial fit, despite the increased size of the 
errorbars. In the lower right panel we show the mean of the fit for the PNGB model (solid line). We recognize that 
the fit is not appropriate for the oscillating PNGB model (dotted line). This behavior does not improve much if 
we apply a second order fit. The polynomial fit hardly can distinguish the reconstructed of the inverse tracker, the 
periodic and the SUGRA potential from a cosmological constant, although their theoretical is completely different 
from W0 = — 1 as evident from fig. |[ 

We will now examine the question of whether we can reconstruct an evolving with the SNe observations. Since 
the values for the iV = 1 fit where sufficient and the errorbars on this fit are relatively small we will concentrate 
in the following on this linear fit. In order to be able to decide if a model is evolving we perform a change of variable 
to more convenient expansion. We can rewrite eqn. (E4) as 



N 



N 



W^{z) = ^^Wi{l + zY = WiZ^ , 



i=0 



(29) 



4=0 



with 



N 

E 

k=0 



Wk ■ 



(30) 
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FIG. 10. The equation of state factor for the periodic modeL The theoretical value is w^{z) is given by the thin dotted 
line. In the top left panel we plot the mean values for the fitted Wij, with the dotted line for the A'' = fit, the solid line for the 
A'' = 1 and the dashed line the for N = 2 fit. The top right panel is the A'' = fit with the shaded region representing the la 
uncertainty levels. The lower left panel is the same plot for A'' = 1 and the lower right panel for N = 2. 



For the A'^ = 1 fit this leads to wq ^ wq +wi and lii ~ wi. The errors in the new expansion coefficients wt can, again, 
be found by Gaussian error propagation 

We calculated these expansion coefficients and their la errors for all the models. In table || we present the expansion 
coefficients for the dark energy models we have discussed. We note that we obtain evidence at the la level for 
evolution for the SUGRA and periodic potential. The only model which is not reproduced correctly in the context 
of an evolving is the brane-inspired model. This is because equation of state factor for the brane model changes 
only relatively early (z > 0.8), where the data does not have so much weight. Note although we get a consistent 
result for the PNGB model, the fit of is relatively badffig. |ri|, lower right panel) and it has with ~ 24.8 an 
unacceptable luminosity distance fit. We also note in fig. Othat the inverse tracker model is not only evolving, but 
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z z 
FIG. 11. The fitted for difFerent dark energy models. The sohd hnes and the dark shaded region corresponds to the mean 
and la error regions of the pure cosmological constant model. The dashed lines and the light shaded region to the periodic 
potential (top left panel), the inverse tracker potential (top right panel) and the SUGRA potential (lower left panel). In the 
lower right panel we show the theoretical (dotted line) and reconstructed (solid line) of the PNGB model. 

also does not behave like a cosmological constant, since « —0.4, which is consistent with the result in table ||. In 
fig. [ij we show the joint probability contours in the wq — wi plane, for the 68.3% and 99% confidence levels. The 
confidence levels for the joint probability are calculated as the regions with = Xmin + with Ax^ — 2.3 and 
Ax^ — 9.21 respectively, which is valid since the errors are symmetric. The shaded region is the 68.3%, resp. 99%, 
confidence region for ~1 < Wcj, — const < 0. In this plot we see that only the SUGRA (thin solid ellipse) and periodic 
potential (thick short dashed ellipse) can be distinguished from a constant at the la level (left panel). At the 99% 
level (right panel ) it is even harder to extract evolving models and only the periodic potential can be distinguished 
from W0 = const. However, if we are just interested in whether a model is evolving we have to concentrate on the 
parameter wi and marginalize over wq, which corresponds to the projection of the confidence region with Ax^ = 1 
for the la errorbars | |89| , |90t . But even if we consider the marginalized errors on wi in table || the only model for 
which we can find evidence for evolution with 99% confidence is the periodic potential. Note that we have omitted 
the PNGB model in this discussion because the reproduced for this model is not valid. 
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wo 


Swo 




Swi 


theoretical evolution 


evolution reconstructed 


A 


-1.00 


0.035 


-0.011 


0.16 






trapped minimum 


-0.99 


0.035 


-0.0057 


0.16 






Brane 


-0.97 


0.034 


0.028 


0.16 


+ 




two exp. 


-0.95 


0.034 


-0.016 


0.16 






periodic 


-0.30 


0.027 


-0.60 


0.11 


+ 


+ 


pure exp. 


-0.84 


0.033 


-0.14 


0.15 








PNGB 


-0.00 


0.025 


-0.94 


0.10 


+ 


+ 


SUGRA 


-0.81 


0.029 


0.31 


0.13 


+ 


+ 


exp. tracker 


-1.00 


0.035 


-0.011 


0.16 






inv. tracker 


-0.40 


0.025 


0.054 


0.10 







TABLE 11. The evolution coefficients with errorbars for the linear fit = wo + Wiz. "+" denotes evolution, "— " no 



evolution and "0" marginal evolution. 
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FIG. 12. The joint confidence regions in the wo — wi plane. In the left panel we show the 68.3% confidence regions and in 
the right panel the 99% region. The shaded region is the uncertainty region for an arbitrary but constant equation of state 
factor with — 1 < w,i, < 0. The thin short dash ellipse is for the trapped minimum models, the thin dot - short dash ellipse is 
for the brane inspired potential, the thin short dash - long dash ellipse for the potential which involves two exponentials, the 
thick short dashed ellipse for the periodic potential, the thick long dashed ellipse for the pure exponential, the thin solid ellipse 
for the Supergravity inspired potential, the thin long dashed ellipse for the exponential tracker solution, and the thick dotted 
ellipse for the inverse tracker. 



Up to now we have used a fixed prior on Clj-n to fit the luminosity distance and reconstruct the equation of state 
factor w^. We will now discuss how our results will change if we have no prior information on Qj^i and use just the 
constraint < Jim < 1. In table III we show the results of the fit with = 1 and no prior on Slni- Note that we 





wo 


5wo 


Wl = Wl 


Swi — 5wi 


f2ni 


SQm 


A 


-0.94 


+ 1.08 

-0.46 


-0.063 


+0.64 
-1.10 


0.31 


+0.07 
-0.31 


inv. tracker 


-0.44 


+ 1.53 

-0.15 


-0.01 


+0.99 
-0.99 


0.38 


+0.23 
-0.38 


periodic 


0.75 


+0.86 
-0.65 


-1.10 


+0.69 
-0.83 


0.41 


+0.08 
-0.21 


SUGRA 


-1.10 


+ 1.31 
-0.20 


0.21 


+0.29 
-1.42 


0.34 


+0.13 
-0.34 



TABLE III. The expansion parameters from eqn. (|24[) for a TV = 1 fit, where Sim is only constrained by < Qm < 1. The 



underlying theoretical models have all f2m = 0.3. 
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show the quantities Wi and not Wi. The reason for this is, that if we include fJm as parameter to fit, the errors on 
the fitted Wi are not symmetric, and therefore not Gaussian. Therefore, we can not perform the error propagation 
using eqn. (|3|). We study two models which are theoretically not evolving, A and the inverse tracker, and two models 
which have an evolving w^, the periodic and SUGRA potential. First, we note that the errorbars on the fitted value 
of r^in are large, and that the mean value, in the case of the periodic and inverse tracker potential is displaced from 
the theoretical value by over 25%! Although, due to the large errorbars the mean value and the theoretical value 
always lie within the la errorbars. In the case of the periodic potential we can recover evolution marginally, since 
wi = — LlOlg gg. However, at the 99% level we can not gain any evidence for evolution. For the SUGRA model we 
can not reconstruct evolution since the obtained value for wi is consistent with no evolution {wi = 0) already on the 
1(7 level. We conclude, therefore, that if we fit simultaneously for ri,„ and wi simultaneously it is poor and has large 
errorbars. 





wo 


5wo 






A 


-0.99 


0.06 


0.30 


0.02 


inv. tracker 


-0.45 


0.11 


0.37 


0.12 


periodic 


-0.28 


0.01 


0.00 (!) 


0.03 


SUGRA 


-0.91 


0.07 


0.38 


0.03 



TABLE IV. The fit results for Q.^ and wo using the Af = fit. 



In table IV we show the results for the constant = fit. First, we note that the errorbars on lup and Jim are much 
smaller than for the linear A'^ = 1 fit. However, we do not expect that this will work for models which are evolving. 
In fact, we recognize that for the periodic potential we get fim = where this is only the "best fit" value because of 
the constraint on Om- If we release the constraint for this fit, Om floats to negative and non-physical values. This is 
because the periodic potential has an evolving Wrj, and an A^ = fit can not reproduce such an evolving model. The 
pure A, the SUGRA and the inverse tracker model seem to give reasonable results, apart from the SUGRA model 
results in a too large value for f2„i, where the true value of ri,„ is outside the la errorbar. This is again due to the 
fact that W(j, is evolving for the SUGRA model. 



In fig. 13 we show the joint probability contours in the wq — Q.^ 



-0.6 




FIG. 13. The joint probabilities for fim and wq for the A'^ = fit for the A model (solid line) and the SUGRA model (dotted 
line). The dark shaded region is the 68.3% confidence level and the brighter shaded region the 99% level. 

plane. The solid line is the pure cosmological constant and the dotted line the SUGRA model. We note that even for 
a marginally evolving model like the SUGRA model that the A^ = fit gives a reasonable result. Note that we did 
not impose a constraint of wo > —1 for the A^ = fit. This constraint basically results in a cut off of the confidence 
regions in fig. ^ below the wq — —1 line. To conclude we remark that we either need a fixed prior on f2ni to be able 
to tell whether the equation of state factor is evolving, or we only consider the constant contribution to which 
would allow us to establish both ilni and accurately. In section VI we will now investigate how different satellite 



17 



designs and priors on fim may change this behavior. 

VI. IMPACT OF EXPERIMENT DESIGN PARAMETERS AND PRIORS ON THE MATTER CONTENT 

We will now discuss the impact of the specific SNAP setup on the estimation on w^. In fig. |lj we show the 
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FIG. 14. The dependence of the statistical errors aij on the SNAP specifications for the linear, = 1, fit. The dotted line is 
(700 ~ Swq the dashed line is — croi = {SwoSwi) and an = Swi the solid line. In the top left panel we show the dependence on 
the number of observed redshifts Nz where the maximal redshift is fixed to Zmax = 2 and the statistical error on the magnitude 
is (Jmag = 0.15 mag. In the top right panel we show the dependence on Umag, where the number of redshifts is fixed to be 
= 2000. In the lower left panel we plot the error matrix as a function of the maximal observed redshift, where we fix the 



number of observed SNe to 2000 and 



0.15 mag. In the lower right panel we show the evolution of the statistical error 



with the maximal observed redshift assuming that 90% of the observations are in the low redshift region < z < 1.2. 

dependence of the error matrix aij on the experimental parameters using the iV = 1 fit for the periodic potential 
model with = 0.3 fixed. The dotted line is Cqo, the dashed line — ctqi and the solid line is CTh. In the top left 
panel we show the dependence on the observed number of redshifts Nz where we keep the magnitude error fixed at 
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fmag = 0.15 mag and the maximal observed redshift is ^max = 2. We recognize that the variances scale as expected 
with 1/Nz, or the errors with l/^/N^. In the top right panel we show that the dependence on the absolute magnitude 
error (Tmag, where the maximal redshift is fixed to Zmax — 2 and the number of observed SNe is — 2000. We note 
that the error matrix is increasing as (Tmag increases as expected. The conservative limit on the total dispersion is 
Cmag = 0.15 mag as stated in table |, which includes using the best currently known methods of standardizing and 
calibrating the luminosity of type la SNe with a residual dispersion of 0.12 mag and a measurement uncertainty of 
0.09 mag after correcting for extinction and using the color of the SNe. Optimistically view that with new methods 
and all the additional information which is available for the SNe with SNAP might reduce the residual dispersion from 
the standardization and calibration to 0.05 mag, although we do not expect to the improvement in the measurement 
uncertainties to be more then 0.08 mag. Therefore, the most optimistic choice results in (Tmag — 0.09 mag. We see in 
the top right panel of fig. |l^ that if we improve the statistical error from (Tmag = 0.15 mag to CTmag — 0.09 mag the 
uncertainties on the fit parameters Wi improve by 70%. This means we are able to reconstruct evolution even for the 
marginally evolving pure exponential model (thick long dashed line, fig. |^) at the 68.3% level with wi — —0.14 and 
Swi = 0.09. 



In the lower left panel of fig. 14 we plot the dependence on the maximal observed redshift where we keep the 
number of observed SNe redshifts fixed to 2000 and (Tmag — 0.15 mag. It is evident from the plot that we can gain 
most accuracy if the observations are done up to a redshift of Zmax ~ 3. Beyond this redshift there is no further 
improvement. In the lower right panel we show more realistic result of increasing the maximal observed redshift on 
the statistical error. In this plot we assume that there is a fixed threshold number An of photons which have to arrive 
in the detector in order to observe an SNe. We neglect the effects of color in this analysis. The brightness or flux of 
a SNe is T = C/Airdf^ and the measured flux in the detector is cx (An) / (AAAt) , where A A is the effective area of 
the detector. From this relation we can work out the time At for which we have to measure the flux from a SNe at 
a particular redshift in order to observe it. We assume that for a fraction p of the total observation time T we find 
SNe in a low redshift region < z < zi and a fraction (1 — p) in the high redshift region zi < z < Zj^ax- In order to 
calibrate the relation we fix the low redshift region with zi = 1.2 and the number of observed SNe in this region to 
Afi = 1850 as in table ^. The number of observed SNe in the high redshift region is then approximately 

jv.^M,'-^ a±itiZ2 

P (l + Zi) + (1 + Z;) (^max - Z;) + (Zmax - 2:/) /3 

In the lower right panel of fig. |l^ we assume that for p equivalent to 90% of the total observation time we discover SNe 
in the low redshift region < z < 1.2. If we choose Zmax = 1.7 we observe Mh = 65 SNe in the region 1.2 < z < 1.7 
in agreement with the values given in table || For redshifts z < 1.2 a similar equation to eqn. (^2|) holds, where in 
this case the entire observation time is spent in the low redshift region. If we observe out to a redshift of Zmax — 2.5 
we can observe Af^ — 48 SNe. From fig. Owe recognize that if we go from Zmax = 1-7 to Zmax = 2.5 we improve the 
statistical error by 10%. We see from table O that such an improvement of the statistical error on wi is only marginal 
and we can not distinguish more models or establish evolution for more models when compared to the case when we 
can just measure out to redshifts of Zmax = 1.7. 

We notice in fig. |ri|, that the "reconstructed" has a region where the errorbars are relatively small, which is 
for the iV = 1 fit in eqn.(p5[) around z/ « 0.2. We find the position z/ of the feature by minimizing 5w^ in cqn.(|2^) 
assuming a symmetric error matrix and a linear fit 

1 + ^/ = -^, (33) 

(Til 



with (Ty the covariance matrix from eqn.(|27j). This feature corresponds to the "sweet spot" in |91| and is also the 
redshift around which pl[ perform the w-expansion in order to obtain uncorrelated expansion coefficients. From this 
expression we already expect that the position of the feature will not vary if we change the number of data points N 
or the statistical error on the magnitude (Tmag and this behavior was established by numerical experiments. In fig. ^ 
we plot the dependence of the position of the feature Zf versus the maximally measured redshift for the periodic 
potential. We plotted the behavior up to the very high redshift z = 20 and recognize that even for such high redshifts 
the feature is still below zj = 0.7. The transition from matter to vacuum domination (when the energy density of 
the dark energy field dominates over the one of matter) for the periodic model is around z w 0.75. So even if one 
could measure SNe at the unrealistic distance of z = 20 we still do not have the smallest errorbars in the interesting 
z = 0.75 region. This behavior is similar for all the dark energy models we studied. 

We examine now the behavior of the statistical error if we fit for the parameters wq and Q,^. In fig. |l^ we show 
the behavior of the covariance matrix aij with a varying Zmax and (Tmag- In this case an is the square of the error 
on the best fit Om- In the left panel we show the evolution with the maximal measured redshift where we construct 
the sampling rate in the high and low redshift bins with eqn. (^) and TVj = 1850, z; = 1.2. We recognize that we 
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FIG. 15. The shift of the minimal error region Zf with the change of the maximal observed redshift z, 
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FIG. 16. Improvement of the statistical error on wo and flm for the fitted A model. The solid line is an = 50,^, the dashed 
line is -~aoi = {SwoSQ,m) and the dotted line is aoo = Swq. In the left panel we show the improvement of the uncertainty on 
the fit parameters due to increasing the maximal observed redshift, where we assume that for 90% of the observation time we 
find low redshift SNe. In the right panel we show the dependence of the statistical error on the residual magnitude dispersion. 



improve the statistical uncertainty by 14% if we measure SNe out to redshifts of Zmax = 2.5 instead of z^ax — 1-7. In 
the right panel we change the residual dispersion (Tmag with the rest of the parameters fixed to the SNAP specification 
in table |. The improvement on the statistical error by going from amag = 0.15 mag to CTmag = 0.09 mag is 64%. We 
will discuss the relevance of this improvement when we examine different priors on the fit parameters wq, wi and fim- 
However we conclude that an improvement of the measurement uncertainties is far more relevant than the ability to 
observe SNe at larger redshifts. 

We will now discuss how different priors on the fit parameter infiucnce the accuracy. We have discussed just 
the most extreme cases of priors on the TV = 1 fit with either Jim — 0.3 fixed or = fixed. In fig. ^ we show 
the decreasing errorbars in the wq — wi plane with different Gaussian priors on fl^. Note that we analyzed the full 
likelihood function and did not assume a Gaussian shape for the probability distribution. The current observations 
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FIG. 17. Error contours of the SUGRA model with different priors on Qm in the wo-wi plane. The solid line contours are 
the 39.3% joint probability regions, which project to la errors on the axis and shaded regions represent the 68% or la joint 
probabilities. The range of increasingly larger contours represents the result using prior knowledge on flm with increasingly 
poorer uncertainty, an,„ = Sflm. The unlabelled, black dotted curve corresponds to the projected-lcr error contour obtained 
for more optimistic dataset specifications and with a prior of 0.25 < Q,m < 0.35. 



provide either just a crude measurement or upper limits on fim. From |5^] we obtain — 0.2t_nj which is too crude 
to result in a significant improvement on our wo-wi-^lm estimation. The X-ray observation in pi| give an upper limit 
of fiin < 0.32 ±0.01 and the Sunyaev Zel'dovich results in are fim < 0-34lg'o3. Future Sunyaev Zel'dovich surveys 
|Eo,^ can possibly determine flm up to an accuracy of 3% if one considers in conjunction with CMB measurements 
p3 |. Future X-ray surveys could determine fim to even higher accuracy independent of any other measurement [p3| . 

We observe in fig. ^ that only very tight bounds of fim ± (^f^m with Jfim < 0.05 improve the statistical errorbars 
which are stated in table If we improve the statistical error on the magnitude to (Jmag = 0.09 and double the 
number observed SNe we can even further improve the accuracy (dotted line in fig. p^ . In table ^ we summarize 



prior (5n„, 


measurement cTmag 






No flrn prior; wi = 


0.15 


0.06 




0.15 


0.15 


0.15 


0.6 


0.05 


0.15 


0.06 


0.2 




0.09 


0.05 


0.12 


0.04 


0.15 


0.05 


0.16 


(fixed Qm) 


0.15 


0.03 


0.12 



TABLE V. Statistical measurement uncertainties on wq and lOi, given supernova magnitude measurement uncertainty, amag, 
and a range of uncertainties, in the independent prior knowledge of (As in fig. ^ , the Supergravity model is used 

here as the example, but the other models give comparable results.) 



our findings for the accuracy of the SNAP measurement with different priors on f2m- If we combine these findings 
with the results in table ^ we see that for the periodic potential the current accuracy on Sim is sufficient to establish 
evolution at least on the la level, however for the SUGRA model we need at least the tight prior with (5ri,„ — ±0.05. 

A further restrictive prior for the linear fit is —1 < w^{z = 0) < which results in the constraint —1<wq-\-wi<0. 
The constraint — 1 < excludes non- minimally coupled scalar-tensor theories [ [76| , p4|j9^ , which we have not included 
in our discussion. If we analyze the results for the = 1 fit with fixed fim this prior does not improve the statistical 
uncertainty. If we do not fix f^m, the constraint —1 < wq + wi < for the linear fit, or — 1 < < for the constant 
fit also does not improve the statistical errors. From this we conclude that either we use a tight measurement of 
ri,n from an independent observation or we can not establish evolution of Wtj, for most of the dark energy models. A 
further prior we have used throughout this paper is the assumption that the universe is flat. For a non-flat universe 
the luminosity distance is given by 
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dh{z) = 



c(l + z) 



S / W + z'f + ^m{l + z'f + + z')3(-o + l)e3-i.' rf^' ^ (34) 



-1/2 



with fifc = 1 — — and = sin(x) for Vtk < and S{x) — sinh(x) for il/j > 0. In fig.ES^ we show the magnitude 




FIG. 18. The difference in magnitudes between non-flat models and the fiducial A cosmology, where the ratio between i},n 
and is fixed. The solid line is for fl^ — —0.2, the dotted line for flk = 0.2, the short dashed line for fl^ — —0.05 and the 
long dashed line for Qk ~ 0.05. The data points correspond to the binned SNAP data as in figj^. 



difference to the fiducial A model for models with non-vanishing curvature ilk, where we fixed the ratio rim/^tj)- The 
solid line is for fifc — —0.2 and the dotted line for 51^ = 0.2. SNAP could clearly distinguish these models on the 
la level. However the current uncertainties on the curvature from CMB and large scale structure observations are 
AQk = ±0.05 on the 2a level |l^j5^. The long and short dashed lines in fig.|l^ correspond to flk = ±0.05 and we 
see that with the SNAP data we can not improve the current constraints. The inclusion of a curvature term in the 
analysis presented in this paper in principle increases the uncertainty of the estimated parameters, however in the 
light that future CMB observations like MAP and Planck will provide even tighter constraints on the curvature (for 
Planck Aflk = ±0.007 |^^), we assume, as mentioned before, a fixed prior of Qk = 0. 

We turn now to the question how systematic errors influence our ability to distinguish dark energy models and 
reconstruct the equation of state factor w^. The systematic error in the luminosity distance is related to that in the 
magnitude CTsys by 



± Sd^ 



(lO±^ - l) 4 



(35) 



Since we marginalize over the magnitude zero point A4 in eqn.(^, we expect the systematic error for the magnitude 
to be zero for low redshifts. We assume a linear drift of the systematic error by 



i^sys 



(36) 



where CTsys is the systematic error at a redshift of z = 1.5. In fig. |l9| we show the influences of the systematic 
uncertainty, where we plot the 68.3% joint probability regions on the estimated parameters with the standard SNAP 
specifications from table ^ The light shaded regions correspond to a systematic error of CTsys = ±0.02 mag, the 
transparent regions to CTsys = ±0.05 mag and the dark shaded region is the result for no systematic error. In the left 
panel we show the result for the TV = 1 fit with a fixed fim. We note that the linear drifting systematic error leads to 
a shift in the wi direction, but the marginalized error on wq changes only slightly. The shift for the periodic potential 
model (dashed lines) is much smaller than for the A model (solid lines). But if we take into account the relative 
sizes of the statistical error depicted by the error ellipses, we find that the shift due the systematic error is roughly in 
agreement with these errors. In both cases the systematic error with CTsys = ±0.05 mag is only in marginal agreement 
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FIG. 19. The effects of tlie linear drifting systematic error on the fitted parameters. The dark shaded region corresponds to 
68.3% joint probability region of the models without systematic error, the light shaded regions to osys = ±0.02 mag and the 
transparent regions to asya = ±0.05 mag. The ellipses with a solid margin are the results for the A model. In the left panel 
we show the results for the = 1 fit with fim fixed. The dashed margins are from the periodic potential model. In the right 
panel the results for the Qn^-wo fit are shown, where the dotted lines correspond to the SUGRA model. 



with the reconstructed values of wi with no systematic error. The value of wi for the A model is lii = —0.011 and 
the one for the periodic potential is wq = —0.60 if we do not include the systematic error. The 68.3% joint probability 
regions for asys = ±0.05 mag only marginally overlap with these values in both cases. However, the 99% confidence 
regions are all in agreement even in the shifted case with the theoretical mean values. We can reconstruct evolution 
of the periodic potential model even for (Tgys = —0.05 mag which results in wi = —0.47 and Swi = 0.11 which is 
still in agreement with wi ^ at the 3a level. We also expect that we can distinguish the same models without 
inclusion of a systematic error since the main difference of the models is along the wq axis as evident from fig. |l^. 
In the right panel we show the results for the ^m-wo fit, where the dotted ellipses correspond to the results from the 
SUGRA model. As before the estimate value of wq is not as strongly affected by the systematic error as the fitted 
value of fiin- Again the 68.3% confidence levels for (jgys = ±0.05 mag (transparent regions) are only marginally in 
agreement with the fitted S7m for no systematic error (dark shaded region). We conclude that the infiuence of the 
systematics is of significant importance for the estimated ilm or wi only at the la level, but is almost negligible for 
the estimation of the constant contribution iuq to the equation of state factor w^. In order to answer the question 
as to whether we can reconstruct evolution and fl^ of the periodic model with the SNAP experiment we have to 
take into account systematic errors. First, we realize that for the standard SNAP configuration from table H with a 
systematic error of (Tgys = 0.02 and a statistical uncertainty we get for the best fit parameter, wq = 0.75lg'g5 ± 0.04, 
wi = —1.07^0 83 i 0.04 and Sim = 0.41^q'2i!'Io o3; where the second error is due to the systematics. We recognize 
that we can only reconstruct evolution, which is the inconsistency of wi with a zero mean, at the la level and that 
the statistical error on is nearly of the order of the mean (50%). The limited amount of linear evolution we can 
reconstruct at the Sa level with this experimental setup is wi « ±2.4 and none of the models we studied has such a 
large linear evolution. If we are able to increase the precision of the SNAP satellite to the limit of a^nag = 0.09 the best 
fit parameters for the periodic model are, wo = 0.75tg;l^ ±0.04, wi = -1.07^5]:^ ± 0.04 and = 0.41+°;°^t°;°2. We 
note that for this setup we can reconstruct evolution at the 2a level, but also that the errorbars on are significantly 
smaller for this setup as for the conservative SNAP specifications from table |. In order to establish evolution on the 
3a level, we need a linear term of the order wi w ±1.5 which none of the studied dark energy models fulfill. In order 
to improve the SNAP proposal even further we assume that we are able to measure twice as many SNe in each redshift 
interval as for the proposed SNAP specifications from table |. For this setup we obtain the following parameters: 
Wo = 0.75lS] jg±0.04, wi = -1.07i[5ji±0.04 and fl^ = 0.4ll^ °^l°;°^. We notice that with the fl^ measurement we 
almost reach the systematic error limit and we obtain relatively tight bounds on the matter contents. Furthermore, 
we have reconstructed evolution for the periodic potential at the 3a level and the limit on evolution we can measure 
with this setup is wi « ±1. So we can conclude that if it is possible to build SNAP with cr,nag = 0.09 and measure 
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around Nz = 3830 SNe out to a redshift Zmax = 1-7 it is possible to establish evolution of the equation of state factor 
Wfj, on the Scr level, if the slope of the linear evolution is \wi\ > 1. 

VII. ALTERNATIVES TO THE SNAP MISSION 

We examine now the question of whether there could be an alternative to SNAP. Therefore, we assume that we 
can improve the low redshift results from SCP and measure about Afi = 160 SNe in the redshift range z — 0.1 — 0.55, 
with a statistical uncertainty of (Tmag = 0.20 mag and a systematic error of agys = 0.05 mag, where the reference 
redshift for the systematic error is z = 0.5. For the high redshift region we assume that we can observe Afh = 100 
SNe in a redshift range z = 2 — 2.5 with the Next Generation Space Telescope (NGST) and the same statistical and 
systematic error as in the low redshift region. In fig. Bin we show the results of these specifications for the A model. 







FIG. 20. The 99% joint probability contours for A'^ = 1 (left panel) and fim-wo fit (right panel), for the SNAP (light shaded) 
and low redshift+NGST (dark shaded) specifications. In the left panel we plot the A model (solid lines) and the periodic 
potential model (dashed lines). In the right panel there is the SUGRA model (dotted lines) and the A model. 



the periodic potential (dashed line) and the SUGRA inspired model (dotted line). In the left panel we plot the 99% 
joint probability regions for the linear fit with a fixed prior on Q^. With NGST specifications we can not establish 
evolution even for the strongly evolving periodic potential model. If we study the constant fit with no constraints on 
Oni (right panel) we can not improve current bounds on H,^ as given in [^9| . We have not included the systematic error 
of this measurement in the analysis, which makes the results even more spread out. Altogether it seems clear that 
a refined low redshift measurement together with the NGST can not provide the same information concerning dark 
energy models as SNAP, although NGST can complement SNAP in the high redshift region. This leaves the question 
open if there are surveys which exploit different physics to gain information about the dark energy content of the 
universe. One possibility are future galaxy cluster surveys | ]93| either with an Sunyaev Zel'dovich or an X-ray survey. 
The analysis for these surveys in the context of dark energy has only been performed for constant models. The 3cr 
errors for the joint probabilities in the Q^-wq plane for a X-ray survey have roughly the same size as for the SNAP 



specification in fig. 13. It is interesting that these methods seem to be complementary to the SNAP observations 
in the sense that the error contours are nearly perpendicular to each other |^^. One of the main drawbacks of this 
method is the unknown evolution of cluster luminosity-temperature relation. How an evolving could infiuence the 
estimation of parameters with galaxy cluster abundance is not clear and should be included in a future analysis. A 
different method was proposed in | ^7t which exploits galaxy counts with the planned Deep Extragalactic Evolutionary 
Probe (DEEP) survey ||9q| . This survey gives roughly the same error contours as the ones we present in fig. |l^. Again 
the analysis of this survey has not incl uded an evolving model. Further possibilities to constrain dark ener gy could 



be the Alcock-Paczynski test [p9,100|, the evolution of density perturbations gravitational waves [101|, lensing 



surveys or observations of the Lyman-a forest |102 
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VIII. CONCLUSIONS 



We have investigated the prospects of future supernova experiments to pin down the nature of the dark energy. We 
have emphasized the importance of choosing a good parameterization scheme for the different dark energy models, 
concluding that parameterizing the equation of state of the dark energy is the most powerful known method, because 
the magnitude-redshift functions m{z) which result give the best fits to most of the dark energy models actually 
proposed. There are two caveats worth noting here. Firstly, there might be an even better parameterization scheme 
out there that has yet to be discovered (that provides even better fits). Secondly, this scheme is unlikely to provide 
good fits to all possible models (the periodic model is an example). The ultimate test for a given model is to simply 
generate m{z) for that model and compare it with the data. 

We then use the proposed SNAP satellite as case study for what might be possible with improved datasets. In three 
recent publications by I. Maor et al. [|6j, P. Astier and our own the prospect of the SNAP mission have been 
briefiy discussed. Maor et al. |8^] argue that SNAP can not distinguish different dark energy models and that it is 
nearly impossible to reconstruct evolution of the equation of state of the dark energy component. Ours and P. Astier's 
p7| results agree with their findings, but as seen in figs, [ij and we show that if we constrain the analysis to a 
constant then it is well within the scope of SNAP to distinguish dark energy models, even if we do not impose any 
priors on fim- If it is possible to exploit the full precision of the SNAP instrument ((Tmag = 0.09 mag) and to constrain 
f2in to 0.05, then it is even possible to reconstruct evolution at the 3ct level as long as the linear evolution today is 
above the \wi\ > 0.6 threshold. These results are confirmed by |^^. As a conclusion we can say that SNAP certainly 
has the ability to distinguish dark energy models from a cosmological constant and possible can put some constraints 
on the evolution of the equation of state of the scalar field component. Whether or not alternative surveys, like a 
X-ray or SZ survey ||9^ or the DEEP survey ||97|] can achieve the same accuracy is currently under investigation. 
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